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[57] ABSTRACT 

A system and method for optimal maximum a posteriori 
(MAP) demodulation using a novel mean phase predictor. 
The mean phase predictor conducts cumulative averaging 
over multiple blocks of phase samples to provide accurate 
prior mean phases, to be input into a MAP phase estimator. 
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MEAN PHASE PREDICTOR FOR MAXIMUM 
A POSTERIORI DEMODULATOR 

STATEMENT REGARDING GOVERNMENTAL 
RIGHTS 

The present invention was made with support from the 
United States Government under Contract Nos. NAS3- 
26959, NAS3-26391 awarded by the National Aeronautics 
and Space Administration (NASA). The Government has 
certain rights in the invention. 

Related Application 

The present application is a continuation-in-part of appli- 
cation U.S. Ser. No. 07/607,582, filed Nov. 1, 1990, now 
U.S. Pat. No. 5,303,269. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to communications systems 
and, more particularly, to phase angle demodulators. 

2. Description of the Related Technology 

Most electronic communication systems in use today 
include a transmitter to transmit an electromagnetic signal 
and a receiver to receive the transmitted signal. The trans- 
mitted signal is typically corrupted by noise and, therefore, 
the receiver must operate with received data that reflects the 
combination of the transmitted signal and noise. Thus, the 
receiver receives data y(t) at a time t, where y(t)=s(t)+n(t), 
the sum of the transmitted signal and additive noise. The 
received data equation can be expanded as follows: 

y(t) = V7 Acos[to 0 i + 0(0] + n(f) 

where A is the signal amplitude, ( 0 o is the carrier or reference 
frequency, 6(t) is the time- varying phase function and n(t) is 
noise. 

Many of these communication systems require that the 
receiver demodulate information in the received signal 
which depends on proper demodulation of the signal phase 
angle at all times during transmission. The demodulation of 
the signal phase angle is problematic in view of the perva- 
siveness of noise. Therefore, for this class of receivers it is 
desirable to optimize phase demodulation, which is equiva- 
lent to optimizing an estimation of the phase function 0(t). 

Analog communication involves modulation that varies 
monotonically with an applied signal, e.g., a change in 
frequency that is proportional to the voltage at the output of 
a microphone as is the case in frequency modulation (FM) 
radio. An analog FM signal can be analyzed by measuring 
the phase of the signal at a sequence of sampling times via 
in-phase and quadrature components, as described in the 
parent patent application, now having U.S. Pat. No. 5,303, 
269, which is hereby incorporated by reference. These 
components can be used to obtain a block of preliminary 
phase samples a(kA), k=l,2, . . . ,N. A corresponding block 
of expected or prior mean phase samples 9 m (kA), k=l,2, . . 
. ,N is also presumably available. The prior mean samples 
are predictions of phase angle based on past data. 

The parent application shows how a Hopfield neural 
network can be used to solve for the maximum a posteriori 
(MAP) phase estimates {0(kA); k=l,2, . . . ,N} based on the 
preliminary phase samples {a(kA); k=l,2, . . . ,N}, the 
corresponding prior mean or predicated phase samples 
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{0 m (kA); k=l,2, . . . ,N}, and the expected autocovariance 
function of the estimated phase shifts, R e [(k-j)A], in a MAP 
demodulator. Ideally, the sample block is at least as long as 
the nonzero support of the autocovariance function. The 
5 most accurate phase estimate is obtained for the sample at 
the middle of the block. If the autocovariance function is 
nonzero over a set of 31 samples between (k-15)A and 
(k+15)A, for example, then a block composed of 31 phase 
shift samples between (k-15)A and (k+15)A would be used 
10 to estimate 0(kA). The output of the Hopfield network also 
contains phase samples between (k-15)A and (k+15)A, but 
samples near the block edges are not as accurate as those 
near the center of the block. Thus, in a MAP demodulator, 
it is desired to increase the accuracy of phase estimates 
15 which are generated based upon blocks of phase samples 
using a MAP estimator and a mean phase estimator. 

SUMMARY OF THE INVENTION 

20 In a maximum a posteriori (MAP) phase angle demodu- 
lator, the present invention utilizes extra MAP phase esti- 
mates for cumulative estimation of the prior mean phase 
vector 0 m . Extra estimates are obtained from overlapping 
data blocks corresponding to a sliding window. Only the 
25 estimate at the center of the window is used as the output 
MAP phase estimate for each window position. 

The presently preferred MAP demodulator, which 
includes a Hopfield network, uses preliminary phase esti- 
3Q mates obtained from data samples {a(kA)=arctan[Q(kA)/ 
I(kA)]; k=l,2, . . . } along with corresponding prior phase 
estimates {0 m (kA); k=l,2, . . . } to compute posterior phase 
estimates {0(kA); k=l,2, . . . } for a block of N samples. 
Prior phase estimates are called “prior mean phase esti- 
35 mates” because they correspond to the mean value of a prior 
probability distribution. The phase values {a(kA), 0 m (kA), 
and 0(kA); k=l,2, . . . } are “unwrapped” by adaptively 
adding 2 ji, 0, or -2 ti so as to reduce large absolute differ- 
ences between succeeding phase values. For example, if 
40 la(kA)-a[(k-l)A]l is greater than n radians, then either 2 tc 
or -2 n is added to a(kA) to produce an unwrapped a time 
series with la(kA)-a[(k-l)A]l§7t. Phase differences 
{a(kA)-9 m (kA)} and {a(kA)-0(kA)} are constrained to lie 
between -7t and n by a similar operation (addition of 2 tc or 
45 —2n so that the differences lie in the prescribed interval). 

The most reliable MAP estimate is at the center of the 
block; the remaining estimated phase samples are extrane- 
ous if the block is translated by one sampling interval A and 
a new set of MAP estimates is obtained for each block 
50 position. The problem is how to exploit the extra, unused 
MAP phase estimates at each block position, as well as how 
to obtain the best possible estimate of {0 m (kA)}. The origi- 
nal patent application suggests using linear (autoregressive) 
prediction of the samples in a given block from MAP phase 
55 estimates in a previous block. 

For highly overlapping blocks, MAP estimates 0(A), 
0(2A), . . . , 0(NA) from the first block position can be used 
as prior mean values for estimates of 0(2A), 0(3A), . . . , 
9(NA) for the next block position, which contains samples 
60 2A, . . . , (N+1)A. Since there is no prior measurement of the 
sample at (N+1)A, the prior value of this sample is predicted 
from 0(A), 0(2A), . . . , 0(NA) via autoregression, where 
0(kA), k=l, . . . , N are determined from the first block 
position. The third block position contains samples 3A, 4A, 
65 . . . , (N+2)A. The prior value of each phase sample 0(3A), 
. . . , 9(NA) can now be predicted by an average of two phase 
estimates obtained from the two preceding block positions, 
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since these samples are common to the two blocks. The prior 
mean value of 0[(N+1)A] is obtained from an average of an 
autoregressive prediction from the first block position and a 
MAP estimate from the second block position. The prior 
mean value of 9[(N+2)A] is obtained via autoregressive 5 
prediction using all the other prior mean values in the block. 
Iteration yields a cumulative average estimate of each prior 
mean phase value in a block. The MAP phase estimate is 
extracted from the center of the block. 

The present MAP demodulator uses a window (block of 10 
samples) with duration equal to the support of the autoco- 
variance function R e [(k-j)A], which is 31 samples in this 
example, with sampling interval A=1 for simplicity. The best 
MAP sample is obtained from the center of the block, but all 
estimated phase samples are used in cumulative averages to 15 
estimate the prior mean phase vector 0 m (k-15), . . . , 
0 m (k+15) for succeeding block positions. 

These and other objects and features of the present 
invention will become more fully apparent from the follow- 
ing description and appended claims taken in conjunction 20 
with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram of one presently preferred 25 
embodiment of a maximum a posteriori (MAP) angle 
demodulator using a Hopfield network. 

FIG. 2 is a functional block diagram of one presently 
preferred embodiment of the mean phase predictor circuit in 
a MAP demodulator, of the present invention, shown in FIG. 30 
1. 

FIG. 3 is a detailed block diagram of one presently 
preferred embodiment of the mean phase predictor circuit, of 
the present invention, shown in FIG. 1. 35 

FIG. 4 is a performance graph plotting root-mean-square 
(RMS) phase error vs. signal-to-noise (SNR) ratio, and 
comparing an arctangent phase estimator, a phase-locked 
loop (PLL) phase estimator, and the maximum a posteriori 
(MAP) phase estimator of the present invention. 40 

DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENTS 

Reference is now made to the drawings wherein like 45 
numerals refer to like parts throughout. 

FIG. 1 illustrates one preferred embodiment of a maxi- 
mum a posteriori (MAP) demodulator 100 . The background 
of a general implementation for the MAP demodulator 100 
can be found in U.S. Pat. No. 5,3103,269. 50 

Essentially, the MAP demodulator 100 receives a phase 
modulated signal such as, for example, an FM broadcast 
signal. A sampling circuit 102 periodically samples the 
integrated products of the received signal y(t) with two local 
oscillator signals, cos(co„t) and sin(co„l), where co„ is the 55 
carrier frequency. The products of y(t)cos(io ( ,t) and 
y(t)sin(( 0 o t) are integrated over an interval A, and the inte- 
grator outputs are sampled at intervals of jA, j=l, . . . ,N, to 
generate sampled in-phase and quadrature components, I(jA) 
and Q(jA), j=l, . . . ,N. The sampled components are fed into 60 
an arctangent circuit 104 and an amplitude circuit 106 . The 
arctangent circuit 104 performs the function arctan(Q/I), 
where the received signal comprises in-phase (I) and quadra- 
ture (Q) components, to provide a time series of phase 
samples a(jA). The amplitude circuit 106 outputs a time 65 
series of sample amplitudes ly(jA)l. The sample phase and 
amplitudes from the circuits 104, 106 are fed into a Hopfield 
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MAP phase estimator 108 . The MAP estimator 108 also 
receives a predicted mean phase 0^(jA). The MAP estimator 
108 provides an estimated phase 0(kA) at its output. Thus, 
the desired goal of the MAP demodulator 100 is to achieve 
an accurate estimate of the phase of the received signal y(t). 

The estimated phase 0(kA) is fed back to a mean phase 
predictor circuit 110, which also receives past phase esti- 
mates. The mean phase predictor circuit 110 generates the 
predicted mean phase, which is fed into the MAP phase 
estimator 108 . 

The presently preferred MAP estimator 108 is structurally 
arranged as follows. A summing amplifier 112 accepts the 
phase samples from the arctangent circuit 104 and subtracts 
each sample from a corresponding predicted mean phase and 
provides the result to a Hopfield network 114 . The Hopfield 
network 114 also receives input signals as a function of the 
amplitude samples after, optionally, processing the signal 
via a scaling circuit 116 , and processing by a shift register 
118 and the multipliers 120 which multiply the signal by the 
autocovariance fijnetion R e [(k-j)A], The output of the 
Hopfield network 114 is fed into a summing amplifier 122 to 
result in the phase estimate 0(kA), the end product of the 
demodulator 100. 

The MAP demodulator 100 processes a block of N 
samples at a time. If the Hopfield implementation converges 
sufficiently quickly (as expected), the block moves by one- 
sample increments. Although the MAP demodulator 100 
produces phase estimates for all samples in the block, the 
sample at the center of the block is the most accurate and is 
passed to the output as the MAP phase estimate. The MAP 
estimate is based on N preliminary phase estimates or phase 
samples {a(nA)=arctan[Q(nA)/I(nA)]}, where {I,Q} are in- 
phase and quadrature components obtained from received 
data. The MAP estimator 108 also utilizes N prior mean 
phase values (0 m (nA)} representing the best available infor- 
mation prior to obtaining each new phase estimate. 

In addition to the preliminary phase estimates {a} and the 
prior mean phase estimates {0 m }, the MAP demodulator 100 
utilizes the phase covariance function R e [(i-j)A], The block 
size N ideally corresponds to the support of this covariance 
function, i.e., the number of samples over which it is not 
identically zero. 

If the Hopfield network has convergence time T C ^A, the 
block will move by only one sample interval A between 
estimates. The phase estimate 9[(k+l)A] will then be based 
on preliminary and predicted phase estimates between 
(k-14)A and (k+16)A. Overlapping blocks imply that nearly 
all the phase samples used to obtain 0(kA) are also used to 
obtain 0[(k+l)A]. Hopfield/MAP phase estimates for 
(k-14)A through (k+15)A are obtained as part of the output 
of the MAP phase estimator 108 for 0(kA). Although the 
estimates near the edges of the block are not as accurate as 
those at the center, all the MAP estimates can be used in a 
cumulative average to obtain an updated version of 0 m [(k- 

14) A] through 0J(k+15)A]. 

The prior mean estimate 0 m [(k+16)A] corresponds to a 
new phase sample at the right-hand side of the block, 
corresponding to translation of the block by one sample. 
This new phase value can be predicted via autoregression 
from the updated values of 0 m [(k— 14)A] through 0 m [(k+ 

15) A] and the autocovariance function R e [(k-j)A]. As the 
sampling window that determines the block of samples 
moves to the right, successive MAP phase estimates at the 
same sample position are obtained and averaged. The prior 
mean phase estimate 0 m [(k-15)A] at the left hand edge of 
the block is based on an average of 30 previous MAP 
estimates and one linear prediction. 
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The predictive process for phase angles, of the present 
invention, is illustrated in FIG. 2. Given a[(k-15)A], . . . 
,a[(k+15)A] and 0 m [(k-15)A], . . . ,0J(k+15)A], the MAP 
estimator 108 yields a set of phase estimates 0[(k-15)A], . . 

. ,0[(k+15)A], with the best estimate 0(kA) at the center. The 5 
new set of phase estimates is used to update the 0 m -vector 
by cumulative averaging of new 0-values with old 0 m - values 
corresponding to the same position in time via a cumulative 
average circuit 130. If N samples 0 x (jA), 0 2 (jA), . . . ,0 w (jA) 
have already been averaged to obtain 0 m iV (jA), the next 
sample is included in the cumulative average by computing 
0n,^iO‘AHl/(N+l)]0 JV+1 (jA)+[N/(N+l)]e m , JV (jA). In a lin- 
ear predictor circuit 132, weights Wj, . . . ,w 31 are obtained 
from the covariance matrix corresponding to R e [(k-j)] as in 
the parent application. 

The MAP demodulator 100, of the present invention, uses 
a window (block of samples) with duration equal to the 
support of the autocovariance function R 0 [(k-j)A], which is 
31 samples in this example, with sampling interval A=1 for 2 o 
simplicity. The best MAP estimate is obtained from the 
center of the window, but all estimated phase samples are 
used in cumulative averages to estimate the prior mean 
phase vector 0 m (k-15), . . . ,0,„(k+15) for succeeding 
window positions. 2 s 

Let 

0 mi (nA)=the update of Q m (nA) corresponding to a block containing 
samples n=j+ 1, . . . , j+N; 


30 


( 1 ) 


The initial (j=0) block contains samples {0(nA); n=l, . . 

. ,N}. For blind initiation of the update process, 0 mO (nA) can 
be set to a(nA), n=l, . . . ,N. For a given transmitted 35 
initiation signal with phase samples 0„(A), . . . , 0 o (NA), the 
initial prior mean phase sample 0 mO (nA) equals 0 o (nA), n=l, 

. . . ,N. 

Map estimates for the first block position {(^(nA); n=l, . 

. . ,N} are obtained using {0 mO (nA); n=l, . . . ,N} and {a(nA); 40 
n=l, . . . ,N}. Updated prior mean phase estimates are then 
computed for samples relevant to the second block position, 
containing samples 2, . . . ,N+1: 


(2) 


45 


The second N-sample block should also include 0 ml [(N+ 
1)A], which can be predicted from the previous N 1 samples 
via the autoregressive estimate 


ff 

0 m i[(iV+ 1)A] = Z fl;9 m i[(iV+ 1 -jW 
]= 1 


(3) SO 


where 


a=[a,a 2 . . . a w l r =« 0 V B 


(4) 


55 


and the elements of the N'xN' matrix R 0 and the column 
vector r e are 


l^ely— ^eE( I- y)A]=£{0[(At+l-OA]} 
[r e ],=i? e (iA)=£{0[(iV+l)A]0[(Ar+l-OA]}. 


(5) 

(6) 


60 


If the total support of the symmetric covariance function 
is N samples, then only 


AT = (N- l)/2, for iV odd 

= Nil, for N even 


(7 ) 65 


previous phase samples have nonzero correlation with 
0[(N+1)A], 

After MAP estimates for the second block position, 
{0 2 (nA); n=2, . . . ,N+1}, are computed from {0 ml (nA), 
a(nA); n=2, . . . ,N+1}, updated prior mean phase estimates 
are obtained for samples relevant to the third block position, 
which contains samples 3, . . . ,N+2: 


0 m2 (nAM 1 /3)02(nA)+e/3)0 ml ( n A), n= 3, . . . ,N 

0 m2 [(A?+l)A]=C/2)0 2 [(W+l)A]+(tA)0 ml [(iV+l)A] 


8m2l(N + 2)A] 


= 2 afimlUN + 2 -J)A] 
T= 1 


(8) 

(9) 

(10) 


where ais defined in a manner analogous to (4)— (6). 

In general, the first prior mean value for sample N+k is 
computed from linear prediction, using samples from the kth 
block position, which contains samples k, . . . ,k+N-l ; 


N 

QmkKN + k) A]= 2 fl/0mi[(Af + k -f)AJ. 
J= 1 


( 11 ) 


For the following block positions, the value of 0„,[(N+ 
k)A] is updated by a weighted sum of MAP estimates and 
past prior mean estimates, until sample (N+k)A is no longer 
part of the block: 

e m ^ +J I(iv+A)A]=[ i/o+i )]e* +x f(A'+A)A]+[y/(,- + 1 [(Aq-*)Aj, 

7=1,2, . . . ,N- 1. (12) 

When j=N— 1, the left hand edge of the block coincides 
with (N+k)A. When j=N, the block no longer contains 
sample 0[(N+k)A], and 0 m [(N+k)A] is discarded by an 
N-sample shift register containing the relevant data block for 
MAP phase estimation. 

The presently preferred mean phase predictor circuit 110 
of the present invention is shown in more detail in FIG. 3, 
in conjunction with the MAP estimator 108 and associated 
interface circuitry. The MAP estimator 108 receives phase 
autocovariance function samples 140, and blocks of phase 
samples and prior mean phases via shift registers 142, 134, 
respectively. The output of the MAP estimator 108 is fed 
into a shift register 146. 

As discussed above, the predictor circuit 110 must be 
initialized with a prior mean phase based upon a linear 
prediction. The samples 140 are also received by the pre- 
dictor circuit 110 via a linear prediction coefficients circuit 
148. A block of initial phase samples 150 is fed into a first 
mean phase circuit 152. The source of the initial prior mean 
phase estimates is described just after Equation (1) above. 
The circuit 152 also accepts signals indicative of prior mean 
phases from the shift register 134. The prediction of 0 m [(N+ 
k)A], performed by the circuit 152, which is the estimate of 
0 m [(N+k)A] for the sample block 0[(k+l)A], . . . ,0[(N+k)A], 
is obtained as in Equation (11) above. 

The circuitry 130 to accomplish the cumulative averaging 
is shown in FIG. 3 as the following components: increment 
and load control 154, prior mean phase register 156, mul- 
tiplier 158, increment and load control 160, phase estimate 
register 162, multiplier 164, summing amplifier 166 and 
current mean phase register 168. The cumulative averaging 
process is described by Equation (12) above for sample 
0J(N+k)A], where Kgl. 

The mechanism for dropping off (discarding) an irrelevant 
prior mean phase estimate at the left edge of a block and 
acquiring (adjoining) a new prior mean phase estimate at the 
right edge of the block is to pass updated prior mean phase 
samples along the length-N shift register 134 corresponding 
to the block of samples that are being processed by the MAP 
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demodulator 108 . When a sample reaches the left edge of the 
block, it is discarded on the next iteration. 

It will be recognized that the present invention can be 
implemented in a variety of different ways. As but one 
example, the demodulator could be implemented using 5 
integrated circuit components, including a significant 
amount of digital logic, but the Hopfield network would 
preferably be constructed using analog circuitry, requiring 
appropriate D/A and A/D conversion, to achieve a fast rate 
of convergence. Of course, standard engineering techniques 1Q 
allow for implementations tending to emphasize computer, 
optical, electrical, or other methods, with engineering and 
economic trade-offs restricting each implementation as a 
practical matter. 

FIG. 4 is a graphical comparison of three demodulators 15 
for angle modulated analog communication signals as mea- 
sured by root-mean-square (RMS) error vs. signal-to-noise 
ratio (SNR). The performance of the MAP phase estimator 
108 in FIG. 2 relative to that of a digital phase-locked loop 
(PLL) and the preliminary arctangent phase estimates a(kA) 2Q 
is shown in FIG. 4. The loop filter in the digital PLL uses the 
autocovariance function R e [(k-j)A] to linearly predict the 
next phase sample on the basis of 31 past phase estimates, 
as in FIG. 2. The predicted phase sample 0 m (kA) is used to 
form a sinusoidal reference signal which is correlated with 
received data over a sampling interval A to produce a phase 
error term 0(kA)-0 m (kA). The PLL estimate 0(kA) is 
obtained by adding 0 m (kA) to the phase error. In an analog 
PLL, this addition operation is performed by an integrator, 
which may be part of the loop filter. 30 

Performance in FIG. 4 is measured in terms of root-mean- 
square (RMS) phase error between a sequence of estimated 
phase samples and phase samples extracted from a simulated 
FM signal. White, Gaussian noise is added to the simulated 
FM signal, and rms error is shown as a function of input 35 
signal-to-noise ratio (SNR). At a given SNR value, the MAP 
phase estimator in FIG. 2 has roughly the same performance 
as a phase-locked loop operating on a signal with twice the 
given SNR. 

While the above detailed description has shown, 4Q 
described and pointed out the fundamental novel features of 
the invention as applied to various embodiments, it will be 
understood that various omissions and substitutions and 
changes in the form and details of the device illustrated may 
be made by those skilled in the art, without departing from 45 
the spirit of the invention. 

What is claimed is: 

1. A maximum a posteriori demodulator, comprising: 
a maximum a posteriori (MAP) estimator receiving a 
block of N phase samples {a(A), . . . ,a(NA)} and a 50 
block of N prior mean phase estimates {0 m (A), . . . 
,0 m (NA)}, wherein the estimator generates a block of N 
phase estimates {0(A), . . . ,0(NA)}; and 
a mean phase predictor, comprising: 

a cumulative averager receiving the block of phase 55 
estimates {0(A), . . . ,0(NA)} and averaging the block 
with corresponding elements from previous block 
averages so as to provide the block of prior mean 
phase estimates {0 m (A), . . . ,0 m (NA)}; 
a linear predictor receiving the prior mean phase esti- 60 
mates {0 m (A), . . . ,0 m (NA)} and predicting a new 
prior mean phase estimate 0 m [(N+l)A]; and 
a shift register receiving the prior mean estimates 
{0 m (A), . . . ,0 m (NA)} and new prior mean phase 
estimate 0 m [(N+l)A] such that the first in time prior 65 
mean estimate 0 m (A) is shifted out and a new block 
of prior mean phase estimates 0 m (2A), . . . ,0 m [(N+ 
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1)A] is provided to the MAP estimator at time 
(N+1)A. 

2. The demodulator of claim 1, additionally comprising an 
arctangent circuit providing the block of N phase samples to 
the MAP estimator. 

3. The demodulator of claim 1, additionally comprising an 
amplitude circuit providing a plurality of amplitude samples 
to the MAP estimator. 

4. The demodulator of claim 1, wherein the MAP esti- 
mator includes: 

a first summer receiving the phase samples and the prior 
mean phase estimates so as to provide a block of N 
phase differences; and 

a Hopfield network receiving the phase differences so as 
to generate a phase error at any one selected sample 
time A. 

5. The demodulator of claim 4, wherein the MAP esti- 
mator additionally includes a second summer receiving the 
phase error and one of the phase samples at the any one 
selected time so as to generate one of the phase estimates at 
the any one selected time. 

6. The demodulator of claim 1, wherein the linear pre- 
dictor comprises a circuit to calculate the equation: 

N' 

e„*[(W + t)A]= X afi^N + k-m 
7=1 

where 

0 mlt is the new prior mean phase block calculated from the 
phase samples beginning at time kA, 

N'=(N-l)/2, for N odd =N/2, for N even, 

a .j is a linear prediction coefficient for sample j. 

7. The demodulator of claim 1, wherein the cumulative 
average comprises a circuit to calculate the equation: 

6 m ,*+/ (N+k)A]=[ l/(/+l)]6 H ,[(Aff*)A]+ [//(/+ 1 )] 1 [ (At+fc) A] , 

1=1,2, . . . ,iV-l 

where 

0 m t k+j is the new prior mean phase block calculated from 
the phase samples beginning at the time (k+j+l)A, and 

Q k+J is the block of MAP phase estimates beginning at 
time (k+j) A. 

8. A method of maximum a posteriori demodulation, the 
method comprising the steps of: 

receiving a block of N phase samples {a(A), . . . ,a(NA)} 
and a block of N prior mean phase estimates {0 m (A), 
• • • ,0 m (NA)}; 

estimating a block of N phase estimates {0(A), . . . 
,0(NA)} from the received block of phase samples and 
the received block of prior mean phase estimates; 

averaging the block of phase estimates {0(A), . . . ,0(NA)} 
with corresponding elements from previous block aver- 
ages so as to provide a cumulative average block of 
prior mean phase estimates {0 m (A), . . . ,0 m (NA)}; 

predicting a new prior mean phase estimate 0 m [(N+l)A] 
from the cumulative average block of prior mean phase 
estimates; and 

generating a new block of prior mean phase estimates 
{0 m (A), . . . ,0J(N+1)A]} at time (N+1)A from the 
cumulative average block of prior mean phase esti- 
mates and the new prior mean phase estimate. 

9 . The method defined in claim 8 , additionally comprising 
the step of repeating the method by substituting a new block 
of N phase samples {a(A), . . . ,a((N+l)A)} and the new 
block of prior mean phase estimates in the receiving step. 
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10. The method defined in claim 8, wherein the generating 
step includes the steps of: 

discarding the first in time prior mean estimate 0 m (A) 
from the cumulative average block; and 
adjoining the predicted new prior mean phase estimate to 
the cumulative average block. 

11. A maximum a posteriori demodulator, comprising: 
means for estimating a block of phase samples and a block 

of prior mean phase estimates, said estimating means 
providing a block of phase estimates; 
means, responsive to said estimating means, for averaging 
the block of phase estimates with corresponding ele- 
ments from previous block averages, said averaging 
means providing a block of prior mean phase estimates; 
means, responsive to said averaging means, predicting a 
new prior mean phase estimate from the block of prior 
mean phase estimates; and 

a circuit receiving the block of prior mean phase estimates 
and new prior mean phase estimate, said circuit capable 
of discarding the first in time prior mean phase estimate 
and adjoining the new prior mean phase estimate to the 
block of prior mean phase estimates. 

12. The demodulator of claim 11, additionally comprising 
a phase sampler providing the block of phase samples. 

13. The demodulator of claim 11, additionally comprising 
an amplitude sampler providing a block of amplitude 
samples. 

14. The demodulator of claim 11, wherein said estimating 
means includes a Hopfield network. 

15. A maximum a posteriori (MAP) demodulator, com- 
prising: 

a mean phase predictor generating a block of prior mean 
phases; 

a MAP estimator receiving phase autocovariance function 
samples R e (iA), a block of phase samples {a[(k+j+ 
1)A], . . . ,a[(k+j+N)A]}, and the block of prior mean 
phases {0 m [(k+j+l)A], . . . ,0 m [(k+j+N)A]}, wherein the 
estimator generates a block of phase estimates {0[(k+ 
j+l)A], . . . 0[(k+j+N)A]}; and 
a shift circuit receiving the block of phase estimates 
{0[(k+j+l)A], . . . 0[(k+j+N)A]} and generating 
updated mean phase estimates that are provided to the 
mean phase predictor. 

16. The demodulator of claim 15, wherein the mean phase 
predictor comprises: 
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a linear prediction coefficients circuit receiving the phase 
autocovariance function samples R e (iA), wherein the 
circuit calculates the linear prediction coefficient a 
using the equation a=R e -1 r e ; 
a first mean phase circuit receiving a block of initial phase 
samples, the linear prediction coefficient a, and a block 
of prior mean phases {0 m [(k+j+l)A], . . . ,0 m [(k+j+ 
N)A]}, wherein the first mean phase circuit calculates a 
to first prior mean phase using the equation 

0 M [(lV+i)Al= Z afinkKN+k-jm 
7= 1 

15 a cumulative averager, comprising: 

a first increment and load control, connected to the first 
mean phase circuit, receiving the first prior mean 
phase 6 m fc[(N+k)A]; 

20 a second increment and load control, connected to a 
current mean phase register, receiving the current 
mean phase 0 m fc+y [(K+N)A]; 
a prior mean phase register receiving the first prior 
mean phase at time a from the first increment and 
25 load control and the current mean phase at time 2A 

from the second increment and load controller; 
a first multiplier, connected to the prior mean phase 
register, receiving prior mean phase 0 m J:+/ _i[(k+ 
N)A} and multiplying by 
30 

i . 

7 + 1 ’ 

a second multiplier, connected to a phase estimate 
35 register, receiving the phase estimates 0* + /k+N) and 

multiplying by 

l 

7 + 1 ’ 

40 

a summing amplifier totalling the results of the first 
multiplier and the second multiplier; and 
the current mean phase register receiving the current 
mean phase 0 m >fct ,[(k+N)A] from the summing 
45 amplifier. 
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